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Galactic cosmic ray (CRs) sources, classically proposed to be Supernova Remnants (SNRs), must 
meet the energetic particle content required by direct measurements of high energy CRs. Indirect 
gamma-ray measurements of SNRs with the Fermi Large Area Telescope (LAT) have now shown 
directly that at least three SNRs accelerate protons. With the first Fermi LAT SNR Catalog, we 
have systematically characterized the GeV gamma-rays emitted by 279 SNRs known primarily 
from radio surveys. We present these sources in a multi wavelength context, including studies of 
correlations between GeV and radio size, fiux, and index, TeV index, and age and environment 
tracers, in order to better understand effects of evolution and environment on the GeV emission. 
We show that previously sufficient models of SNRs’ GeV emission no longer adequately describe 
the data. To address the question of CR origins, we also examine the SNRs’ maximal CR contri¬ 
bution assuming the GeV emission arises solely from proton interactions. Improved breadth and 
quality of multiwavelength data, including distances and local densities, and more, higher resolu¬ 
tion gamma-ray data with correspondingly improved Galactic diffuse models will strengthen this 
constraint. 
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1. Introduction 

In this work we will shortly describe the results obtained from the systematic analysis of 
supernova remnants (SNR) in 3 years of F<?rm/-Large Area Telescope (LAT) data. We will correlate 
our spectral and spatial results with those obtained in other wavelength, trying to better understand 
SNRs’ emitting mechanisms, ages and environments. At the end of this work we will discuss under 
some assumptions the contribution of the SNRs to the observed cosmic rays (CR). The method of 
this analysis is described in [ 1 ] while the full analysis will be shown in a future publication [ 2 ]. 

2. Summary of results 

Using three years of data and starting from the spatial information of 279 SNRs [3] we de¬ 
tected 102 candidates with a final source TS > 25. Of the 102 detected candidates, 36 passed 
our association probability threshold. Of these, 30 SNRs 11% of the total) show significant 
emission for all alternative interstellar emission models (IBM, see also [4]) and are classified as 
likely GeV SNRs. An additional four were identified as sources which are not SNRs; two other 
candidates were demoted to marginal due to their dependence on the IBM. Of the sources likely to 
be GeV SNRs, 17 show evidence for extension (TS^^r > 16). Only sources associated with SNRs 
G34.7—0.4 and G189.1+3.0 show evidence of significant spectral curvature in the 1 — 100 GeV 
range and are fit with logParabola spectra. Of the classified candidates, four extended and 10 point 
SNRs are new. Bor all the sources we evaluated their spectral and spatial characteristics with sys¬ 
tematic and statistical uncertainties. We reported 245 Bayesian [5] upper limits (UL) for all the 
non-detected SNRs. In Bigure 1 the observed spectral parameters of all the SNRs are shown: the 
photon flux range is larger than two order of magnitude while the power-law (PL) index range from 
1.5 to ~ 4.0. 

3. Multiwaveleugth correlatious 

We find that the best-fit GeV diameter is within errors of the radio diameter for most of the 
candidates classified as associated with an SNR, as shown in Bigure 2. This result support the 
hypotheses which have both GeV and radio emission arising from the same location, e.g. [ 6 ]. 

Bigure 3a shows the flux from synchrotron radio emission at 1 GHz in comparison to the 7 - 
ray flux at 1 GeV. A correlation may exist between the radio and GeV flux [ 6 ], if the same lepton 
population is directly producing both emissions or if they scale with some underlying physical 
parameter such as ambient density. In [2] we do not find evidence for a significant correlation, but 
we also cannot strictly rule out an intrinsic correlation for various effects (e.g. small energy range 
used for the GeV analysis, possible change in the spectral index etc). 

Bigure 3b compares the deduced radio spectral index a with the 1 — 100 GeV photon index 
r, theoretical predictions for two emitting mechanism are shown on the plot (see label). Nearly 
all candidates have softer-than-predicted 7 -ray photon indices given their radio spectra. The three 
young SNRs in blue are most consistent with a single underlying particle population, and it has 
been suggested they emit via inverse compton (IC, dashed line) at GeV energies. The observed 
soft GeV spectra relative to the radio has several potential explanations (e.g. different PL index for 
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Figure 1: The distribution of fitted photon index and flux in the energy range 1 — 100 GeV. Open 
circles indicate extended SNRs while filled circles indicate point-like sources. All SNRs that passed 
classification are shown as black unless also classified as young non-thermal X-ray SNRs (blue) 
or as interacting with MCs (red). Candidates which did not pass classification but still had both 
fractional overlaps > 0.1 are grey. If they are also young or interacting, they are outlined in blue 
or red, respectively. Statistical error bars have caps; error bars without caps present the systematic 
error, described in [1]. The same color code is used in all the following plots. 


leptonic and hadronic populations, the emitting particle spetra has breaks or different zones with 
different properties dominating the emission at different wavelengths). 

In Figure 4 we plot the PL index in the GeV versus TeV range for all SNRs observed, the ma¬ 
jority of them seems to have TeV indices that are softer than their GeV indices. This characteristic 
suggests spectral curvature, potentially reflecting a change in spectral slope of the underlying parti¬ 
cle population(s’) index or indices, or that another particle population dominates at higher energies 
or the emssion mecchanism is different. 

4. Age Vs. Environment 

In Figure 5a the SNR ages from literature versus the 1 — 100 GeV photon indexes are shown. 
Young SNRs tend to have harder GeV photon indices than interacting SNRs, which are likely 
middle aged, though the scatter in age for the two classes is one to two orders of magnitude. This 
trend may be due to the decrease of the maximum acceleration energy as SNRs age and their shock 
speeds slow down or the interacting SNRs may be more luminous due to their interactions with 
denser surroundings not yet reached by younger SNRs. Ultimately, the large scatter observed in 
luminosity will likely reflect effects due to both age and environment. 

To investigate the role of environment in the trends for the young and interacting SNRs, we 
examined the GeV luminosity versus radio diameter in Figure 5b. The square of the physical 
diameter (D) can be regarded as a reasonable indicator for SNR age and environment, see [2]. We 
note that any apparent correlation between the luminosity and may be due to their inherent 
dependence on distance (squared). Figure 5b shows that, for the detected candidates, interacting 
SNRs are generally more luminous for a given physical diameter than young SNRs, though there 
is large scatter. This suggests that SNRs at the same age may be more luminous because they have 
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Figure 2: The radio diameters of the SNRs from Green’s catalog are correlated with the fitted GeV 
diameters for those candidates with significant extension. The solid line represents equal radio and 
GeV diameters. All cases of detected extension have diameters greater than 0.2°. The ticks denote 
the radio extension of GeV point-like candidates, colored in order of their characteristics (young 
or interacting) and by their classifications (well defined or marginal). The small ‘x’s bracketing 
the points show the minimum and maximum GeV extensions allowed such that the source remains 
classified or marginally classified given the radio position and extension and best fit GeV position, 
see [2]. 

encountered denser gas. With the addition of upper limits, we find some interacting candidates are 
constrained to lie below the luminosities of most young SNRs. Thus, as we continue to detect SNRs 
with increasingly fainter y-ray fluxes, we are likely to find less separation between the luminosities 
of the two classes. 

5. Constraining Cosmic Ray acceleration 

This systematic search for y-ray emission from known SNRs provides an opportunity to un¬ 
derstand if Galactic SNR population is able to accelerate and release CRs with the appropriate 
composition and flux up to the transition between the Galactic and extragalactic components. In 
the following, we assume that the y-ray emission from SNRs probed with LAT entirely arises from 
the interaction of CR protons and nuclei with the surrounding ISM or circumstellar medium through 
the production and subsequent decay of Given that two other emission mechanisms involving 
accelerated leptons, namely non-thermal bremsstrahlung and IC scattering, could also contribute 
in the y-ray domain, the constraints derived from the LAT measurements should be considered as 
upper limits on the CR energy content in SNRs. 

Figure 6a shows the y-ray flux from p-p interactions [7] in the (£’cR,max, Tcr) plane and shows 
its dependence on F^cR,max for different CR spectral indices, with J = 1 kpc, n = 1 cm“^, Eqr = 
0.01, and F^sn = 10^^ erg (for further details see [2]). The CR energy content is computed for 
particle momenta above 10 MeV and is related to the total supernova explosion energy with the 
relation £’cr = ^cr^^sn- The y-ray flux is nearly independent of the CR maximal energy as long 
as £’cR,max ^ 200 GeV and Fcr ^ 2. In this case, the y-ray flux can conveniently be approximated 
and parametrized as shown in Figure 6b. Details on this parametrization and references can be 
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(a) Flux comparison. (b) Index comparison. 


Figure 3: Left: comparison of y-ray and radio spectral flux densities for all SNRs and candidates. 
For all SNRs that were not detected or which failed classification, grey triangles indicate upper 
limits at 99% confidence, computed assuming the radio location and extension. Right: comparison 
of radio spectal index, a, and GeV photon index, F. The expected correlations for decay or 
e^ bremsstrahlung (solid) and inverse Compton (dashed) are plotted. Ticks along the right hand 
side of the right plot show the 1 — 100 GeV photon indices of those SNRs without reported radio 
spectral indices. 


found in [2]. We use the relationship between an SNR’s y-ray flux, density, and distance shown in 
Figure 6b to determine the maximal CR energy content Fcr through 8 cr contributed by every SNR 
for which we have measured the y-ray flux and photon index or derived an upper limit at the 95% 
confidence level. In the case of a detected SNR, the photon index, known to reproduce the spectral 
shape of the parent CR proton/nuclei spectrum above a photon energy of 1 GeV [7], is taken to be 
equal to Fcr. In the case of the upper limits, we assume an index of 2.5 (i.e. the average value of 
the detected SNRs). We use the canonical value of 10^^ erg for Fsn- 

Translating SNRs’ y-ray measurements into constraints on their contribution to Fcr also re¬ 
quires knowledge of their distances and effective densities, both quantities are taken from the lit¬ 
erature. Figure 7 shows the constraints on the CR energy content for the population of known 
Galactic SNRs. As is clearly visible for the first two subclasses of SNRs, the estimates and upper 
limits on the CR energy content span more than three orders of magnitude, from a few x 10"^^ erg to 
several x 10^^ erg. For the interacting SNRs that lie above the 8 cr = 1 dashed line, the densities ex¬ 
perienced by the CR particles in the molecular cloud (MC) interaction region are likely much larger 
than those used in this analsysis. Interacting candidates’ lying above this limit suggests that they are 
likely the sites of hadronic interactions in dense environments. In contrast, most of the young SNRs 
lie at or below this luminosity limit, suggesting that IC processes may contribute to their measured 
luminosity. New information about distances and densities can give us further insight on the ability 
of known SNRs to provide the observed CRs. The usual assumption of 8 cr = 0.1, required in order 
for the Galactic SNR population to supply the CR flux observed at Earth, is compatible with the 
results of the Fermi LATSNR catalog, for more details see [2]. 
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Figure 4: GeV index compared to published index measurements from Imaging Atmospheric 
Cherenkov Telescopes (lACT). The line corresponds to equal index values. The ticks represent 
the GeV candidates with indices in the range of those with a TeV counterpart but with no TeV 
measurements themselves. Yet many of these are hard: 12 candidates and 10 marginal candidates 
have indices harder than 2.5, suggesting they may well be observable a TeV energies. 




(a) Age versus GeV spectral index (b) GeV luminosity versus (radio diameters)^. 

Figure 5: Left: Age versus GeV spectral index. For those with ages in the literature, the young 
(blue) SNR candidates are separated in this phase space from the identified interacting candidates 
(red). The ticks on the right show indices for GeV candidates without well-established ages. Right: 
The 1 — 100 GeV luminosity is plotted against the square of the radio diameters in pc of those SNRs 
with known distances. 


6. Conclusion 

We examined our GeV candidate SNR population in light of multiwavelength observations 
in order to better understand both SNRs’ characteristics and potential for accelerating hadrons. 
In this work we have studied correlation of our results with spectral and spatial information from 
radio and TeV bands, in order to understand if the underlying CR population is the same and has the 
same spectral characteristics. Within the limits of existing MW data, our observations generally 
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Figure 6: Left: Under standard assumptions (see text), an SNR’s y-ray flux in the 1 —100 GeV range 
can be related to the accelerated CRs’ maximal energy £'cR,max and spectral index Fcr for a given 
CR energy content above a particle momenta of lOMeV (Pcr = £^cr/^sn = 0.01), effective 
density (1 cm“^), and distance to the SNR (1 kpc). Right: the relationship between the SNR’s y-ray 
flux in the 1 — 100 GeV band and Fcr, max for different values of Fcr. For Fcr, max ^ 200 GeV and 
Tcr ^ 2, the flux is weakly dependent on the CR maximal energy. 



Figure 7: Estimates of the CR energy content (in units of 10^^ erg) for all Galactic SNRs, divided 
into three categories each sorted by Galactic longitude: SNRs with existing distance and density 
estimates (upper left panel); SNRs with known distances (upper right panel); and SNRs with un¬ 
known distance and density (lower panel). The two dashed lines indicate a CR energy content of 
10 and 100% of the standard SN explosion energy. 
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support previous findings of changes in spectral slope at or near TeV energies and a softening 
and brightening in the GeV range with age, yet we see indications that new candidates and new 
multiwavelength data may provide evidence of exceptions to this trend. In examining tracers of age 
and environment, we found that the classified candidates followed the previously observed trends 
of young SNRs being harder and fainter at GeV energies than older, often interacting SNRs. 

We find, also, that the limits on CR energy content span more than three decades, including 
many interacting candidates for which the densities in the interaction regions are much greater than 
the nominal density assumed in the calculation, and young candidates with efficiencies below the 
nominal ^ 10%, consistent with possible leptonic emission predictions (e.g. IC). The contribution 
from all SNRs, particularly those with flux upper limits, is beginning to constrain the energy content 
put into CRs from the known SNRs to less than 10%, particularly in regions of well characterized 
IBM background. We find that the candidates and upper limits are generally within expectations if 
SNRs provide the majority of Galactic CRs. 
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